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Abstract. - We analyze the coupling of qubits mediated by a tunable and fast element 
beyond the adiabatic approximation. The nonadiabatic corrections are important and even 
dominant in parts of the relevant parameter range. As an example, we consider the tunable 
capacitive coupling between two charge qubits mediated by a gated Josephson junction, as 
suggested by Averin and Bruder. The nonadiabatic, inductive contribution persists when the 
capacitive coupling is tuned to zero. On the other hand, the total coupling can be turned off 
(in the rotating wave approximation) if the qubits are operated at symmetry points. 



Introduction. - Most approaches to quantum computation rely on tunable interactions 
between pairs of qubits. These should be turned off for independent single-qubit manipulations 
and to reduce crosstalk of qubits during readout, and they should be turned on selectively in a 
controlled way in order to perform two-qubit logic gates. In principle, refocussing techniques, 
known from NMR, can be used to suppress the effect of fixed (uncontrolled) couplings, but 
they require considerable overhead and precision of the pulses. Thus tunable setups, which 
allow minimizing the residual couplings, are desired. At this stage, most experiments with 
Josephson qubits [1-6] have been performed with direct and fixed couplings. But various 
schemes with tunable couplings have been proposed [7-20]. Some of them gain their tun- 
ability from ac-driving [19,20], but the majority relies on additional circuit elements, such as 
switchable Josephson junctions, inductors, LC circuits (cavities), or further qubits. One can 
distinguish two operation principles: For "resonant" couplers [9, 10, 15] the coupling element, 
typically an oscillator, is tuned into resonance with one or both qubits. Alternatively one can 
use "adiabatic" couplers [7,8,11-14,16-18], where the coupling element has a much higher ex- 
citation energy than the qubits and remains in its ground state while mediating the coupling. 
In the following we will concentrate on these adiabatic coupling schemes. 

Some early proposals [7, 8] made use of a fixed coupler, but gained tunability by using 
SQUID-type qubits with flux-controlled Josephson energies. They suffer from the difficulty 
that the coupling can be switched off only if one succeeds in fabricating identical junctions. The 
alternative approaches, proposed more recently, employ a tunable coupler which - ideally - can 
be tuned to cancel fixed existing couplings [11,14,18]. In addition, by modulating the coupling 
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constant around such "zero" points one can tune two qubits with different energy splittings 
into resonance, while keeping both at their symmetry points where decoherence effects are 
minimized [17,18]. In this Letter we will show that nonadiabatic corrections are important 
around such "zero" coupling points. For example, we find that a gated Josephson junction 
produces, in addition to a tunable capacitive coupling, an inductive one, which dominates 
when the capacitive interaction is switched off. The importance of inductive corrections in 
the charging regime of Josephson junctions has recently been pointed out also by Zorin [21]. 
Certain nonadiabatic corrections were also noted in ref. [18]. 

Below, we first present a general theory of tunable adiabatic coupling. As specific example 
we analyze the setup with tunable capacitive coupling proposed by Averin and Bruder [11] 
and demonstrate the importance of inductive corrections in different regimes. When the 
charging energy of the coupling junction dominates over the Josephson energy, Ec ^ Ej, 
the inductive coupling is weak, but is important when the capacitive interaction is tuned to 
zero. In the opposite limit, Ej » Ec, the inductive coupling always dominates over the weak 
capacitive interaction, and we recover earlier results for inductively coupled charge qubits [22]. 
In the present discussion we concentrate on charge qubits, but our general formulation and 
conclusions apply equally to other types of qubits and couplers. 

General theory. - We consider a system of two qubits coupled via an adiabatic coupler 
and described by the Hamiltonian 

H = Ho + V = ifqubits + -ffcouplerd) + V . (1) 

The coupler is controlled by a parameter 7. Wc assume an interaction of the form V = \AB, 
where A is the coupling constant, A an observable of the coupler, and B an arbitrary function 
of observables of both qubits. In general iJqubits includes a direct and fixed coupling between 
qubits. In the following wc consider the situation, where this fixed coupling is (nearly) canceled 
by tuning the coupler appropriately. 

The coupler Hamiltonian can be brought into diagonal form, iJcoupier = J2m=o -^"^ 1"^) ("^1 > 
where both the cigcnstatcs |m) and energies Em depend on 7. The small parameter govern- 
ing the adiabatic approximation is |-ffqubits|/(-E'i — -^o), where |i?qubits| denotes the maximum 
energy difference between the eigenstates of -ff qubits- To proceed we integrate out the cou- 
pler and derive an effective Hamiltonian in the subspacc of slow degrees of freedom of the 
qubits. The effective interaction Vge is found by projecting the full time-evolution operator 
S{t,0) = T exp [ — i Jg Vi{t')dt' /h] (in the interaction representation) onto the ground state of 
the coupler, 5eff(t,0) = (0| S{t,0) |0), and defining 

Sesit, 0) = T exp [ - V,sjit')dt'] . (2) 

We assume the interaction to be weak, such that \{A){B)/{Ei — Eq) ^ 1, for all matrix 
elements of A and B. Expanding the evolution operator Ses{t,0) up to second order we find 

Ses^l - J j\o\Ai{t')\0)Bi{t')dt' 

+ (^j) J^\o\Ai{h)Ai{t2)\0)Bi{h)Bi{t2)dhdt2 . (3) 
The first order contribution to I4ff,i is thus given (in the Schrodinger representation) by 



V^^>=X{0\A\0)B 



(4) 
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To find the second order contribution we insert the unity, |m) (to|, in the last term of 
eq. (O between the two coupler operators. The part with m = is the second order, reducible 
term in the expansion of T exp [ — i /q K:ff\(^')'^^'/^]- The terms with m 7^ give rise to a 
new contribution V^ff j to the effective interaction Hamiltonian, defined by 



\2 °° 



J2\{0\A\m)f-I„ 



(5) 



m— 1 



where /,„ = J* J*' Bi(ti)Bi(t2) e'(^o-^'")(*i-*=^)/'* dti dia- Introducing t = h - t2 and T = 
{ti + ^2)72 we rewrite the integral as 

ft rF(T) 

Im= dT dTBi(r + T/2)Bi(r-T/2)e-'"-''^ (6) 

Jo JQ 

where LUmO — {Em — Eo)/h and F{T) = 2 min(T, t — T). Because of the assumed separation 
of time scales, Bi varies slowly on the time scale w~0' '^^^ ^'^ expanded, 

t .F{T) 

dT / dre-'"-"^ 



B,{T) + -BiiT) + —Bi{T) 



B,{T)-^B,{T) + ^B,{T) 



(7) 



For T ^ i-^mo '-^^^ F{T) 00. Substituting the resulting I„i into eq. (jsj we find the 
second-order contribution to the effective Hamiltonian (in the Schrodinger picture) 



V 



(2) 



A' 



off 



-^|(0|^|m)|^ 

7n—l 



B^ 



i[B,B] 



2^2 - {B,B} 



(8) 



Note that the time derivatives in the Schrodinger picture should be understood as commuta- 
tors, e.g., B = i[77qubits, -B]/^- 



The effective Hamiltonian so far is HcS ~ H, 



term of V„ 



V 



(2) 



qubits 



y{2) 



To simplify the last 



(2) 
iff ■ 



an "integration by parts" can be performed which adds a full time derivative to 
This is achieved by a unitary transformation Heft = UHcffU^^, where U = exp [iQ] and 
I (0 \A\ m) I' {B, B} . This yields V^l^ = V^l^ 



m) 



(2) 

eff 



A2 



{B, B} . This yields V^^' = V^^ 
B^ i[B,B] 



i[Q, iJqubits], i-e. 



n ^ 

m— 1 



|(0|A|m) 



The full effective Hamiltonian thus reads 



i?cff = i/qubits + + ^eff' + O(A^) . 

Together with eqs. and © it constitutes our main result. The tunable part of the coupling, 
eq. ©, depends via |to(7)) and uJmO^ci) on the parameter 7. For time-independent qubit 
operators Bi we recover the results of the Born-Oppenheimer approximation. Indeed, eqs. Q 
and © (with 5 = 0) can be obtained by expanding the ground-state energy Eo{B) of the 
coupler and the interaction term XAB in A. In general Bi(T) depends on time, and we 
obtain nonadiabatic corrections. In the considered limit, they are smaller than the adiabatic 
contributions, however, they still dominate the total qubit coupling, if the tunable adiabatic 
part in eq. ^ (approximately) cancels further direct coupling terms in -ff qubits- Below we 
analyze the nonadiabatic corrections for specific systems. 



.(2) 



B' 



(9) 



(10) 
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Fig. 1 - Tunable capacitive coupler proposed in ref. [11]. We slightly modified the system by in- 
troducing the gate capacitor at the coupler. This influences the expressions for the gate charges 
but does not change the form of the Hamiltonian. The effective qubit Josephson energies can be 
tuned using control fluxes in the SQUID loops. To simplify expressions, we assume the SQUIDs 
to be symmetric. Then Ejk ~ 2i5jj. cos(7r<l>xfc/$o) and C.jfc = 2Cjfc, where Ej^. is the Josephson 
energy of one of the SQUID's junctions and Cj^ is its capacitance. The charging energies are given 
by Ec = 2eV(CE - Efc ClJCsk) and Eck = ^e^Csk, where Cs = Cj + Cg + Cn.i + Cm2 and 

Csk = Cjfe + Cgfe + Cmk ■ 



Tunable capacitive coupling. - The system shown in fig. ^ which is similar to the one 
proposed in ref. [11], provides an example of tunable capacitive coupling. It is described 
by the Hamiltonian H — Ec[n — rig — q{ni, 712)]'^ ~ Ej cos0 + J2k=i 2 with single-qubit 
Hamiltonians Hk = Ecki^k — "-gfc)^ ~ Ejk cos (f>k- Here (f> and are the phase differences 
across the coupling junction and of the qubits, respectively. The (dimensionless) gate charges 
of the qubits only depend on the apphed gate voltages, rigk = CgkVgk/2e. But the gate charge 
of the middle junction consists of rig — CgVg/2e plus the gate charge induced by the two 
qubits, g(ni,7i2) = - X]fc(Cmfc/Csfc)(nfc - Ugk). 

In order to connect to eq. we rewrite the Hamiltonian of the system in the form 

-ffqubits = ^ Ecq^ + 2Ecngq , 

-ffcouplor = Ec{n - Ugf - Ej COS(/) , 

V = -2Ecnq . (11) 

Thus we have A — n, B = q{ni, 712), 7 — rig, and A = —2Ec, and the small parameter of the 
perturbation theory is CkEc / {fu^Jmo) with the small constants Cfc = —Caik/C^k entering the 
operator q{ni,n2). Proceeding as outlined above we find the effective qubit Hamiltonian 

°° 2E'^ 

ffeff = ^ i/fe-2ii;c(0|7i-ng|0) ^ -^|(0|n|m)|2-i[q,g] 



k=l,2 m=l 



mO 



\ m— 1 / m— 1 '^^^ 

where q — — J2j ~K''^3 sin^j and [q, q\ = i^^ cos (f>j. By separating the single-qubit and 

coupling terms we arrive at H^s — X]fc=i 2 ^'k + ^int- The single-qubit terms iJ^ = + ^Hk 
acquire a small correction, bEl^^ — 0{ck), while the interaction has two contributions, 

^^int = Xcing)nin2 + Xiirig) sin^i sin 02 , (13) 
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with 

\ m — 1 ^^^^ / 

E -^|(0|n|m)n ciC2i?ji£;j2 . (14) 

m=l "^"O / 

The first term (Ac) corresponds to a tunable capacitive coupling, while the second one (Aj), 
which arises due to the nonadiabatic corrections, corresponds to an inductive interaction, 
coupling the operators of current flowing into/out of the qubit islands. This inductive coupling 
is weak due to a small factor EjiEj2/ (fi^^mo)'^ ■ Yet, as we show below, it dominates, when 
the capacitive interaction vanishes or is small. 

Coupling junction in the charge regime, Eq 3> -Ej. - In this limit, working close to the 
degeneracy point rig = 1/2, we can approximate the coupler by a two- level system with the 
two charge states \n = l),\n ~ 0) as basis. The coupler Hamiltonian then reads 

-^coupler = --Ec(2ng - 1)0-2 - 7;Ej(Tx = --hioiQPz. (15) 



In the second, diagonal form we used fiwio = \/ + E'^{2ng — 1)^ , sinyy = Ej / {huio), 
and ffz = cosr] pz — sinryp^.. Expressing the non- vanishing matrix element in eq. p4|l as 
(0 \n\ 1) = — ^ sin 77 we obtain the coupling constants 

/ E(j \ 2£'2 
Ac(ng) = 2Ec ( 1 - ^ sin^ V ] C1C2 , Ai(ng) = -— ^ sin^ rj ciC2£'ji£'j2 • (16) 



When the capacitive interaction is switched off, at sin'^ r/(ngo) = Ej /Ec, the inductive coupling 
persists 

Ai(ngo) = -2ciC2{Ec/Ejy/^iEnEj2/E,i) . (17) 

Assuming that the ratio Ec/Ej is not extremely large, and using Ej > Ejk (which is required 
in the charge regime to fulfill the adiabaticity condition), we find |Ai(ngo)| ^ ciC2i?jfc. Further 
away from the limit Ec ^ Ej, i. e. for Eq ^ Ej , the coupling constants can still be calculated 
numerically using eq. (|14|l . An example for Ec/Ej = 2 is shown in figs.j^J^) and b). 

Coupling junction in the phase regime, Ec <^ Ej. - In this limit the phase of the coupler 
(j) remains small, {(jP) ^ 27r. Then we can expand Ejcoscj) « Ej — Ej(jP/2, and the coupler 
Hamiltonian 

i?couplcr=^-Bc(n-71g)2+^jc/.V2 , (18) 

reduces to that of a shifted harmonic oscillator. Making use of n — Ug = i{a'^ — a) hioiQ / {AEc) 
with uoiQ — j^y'2EjEc , we evaluate Ac and Ai from eq. H14(l . It turns out, that in the 
considered limit and approximation the capacitive part of the interaction vanishes while the 
inductive coupling constant, Ai(ng) = — ciC2(i?ji£'j2/£'j), is independent of Ug. A more 
precise analysis of the Hamiltonian of the coupler yields an exponentially weak rig-dependent 
capacitive interaction [11], but the inductive coupling still dominates. In this limit the setup 
reduces to that considered in ref. [22], with the Josephson junction of the coupler playing 
the role of the inductance of the LC circuit with L = ^'^/{ATr'^Ej). Results for the regime 
Ec ^ Ej are illustrated in fig. |21c). 
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Fig. 2 - Coupling constants Ac(ng) and Ai(ng) in units of C\C2Ec obtained from eq. 11411 . a) Ec = 
2Ej, E,n — E]2 ~ 0.2E,T. b) A closer look at the vicinity of the "zero" point, c) Ec — 0.5Ej, 
Eji = Ej2 = 0.2Ej. For other values of E.n and Ej2 the inductive coupling constant can be obtained 
by simple scaling (see eq. Ill4ll 'l. 



Switching off the coupling at the symmetry point. - Eq. (|14|l or eq. H16|l show that 
the coupling cannot be switched off completely. Yet, in certain situations, at least in the 
rotating wave approximation (RWA), the most important part of the coupling can be switched 
off. As an example we consider two qubits (fc = 1,2) in the charging regime with Hk = 
-^BzkC^zk - ^BorkCTxk- Ignoring renormalization effects we have B^k = Ecki^Ugk - '^)crzk 
and B^k = Ejk] renormalization effects introduce a gate charge dependence due to voltage 
crosstalk, S^j, — Bzk{ngi,ng2,ng). The effective interaction Hamiltonian can be written as 

AcK) Ai(ng) 

-nint — -. Crzl<7z2 H -. <^%i\^y2 ■ 



Diagonalization of Hk leads to Hk = —^AEkPzk with AEk — a/B^j, + and azk = 
cosrjkPzk — sinrjkPxk with sinrjk = Bxk/ AEk- The interaction Hamiltonian then reads 

Hint = ^^^^^(cOS77iP2i - sin77iP2,i)(cOS 772Pz2 - Sa^V2Px2) + ^'^^^^ PylPy2 ■ (20) 

If we assume the qubits to have equal or similar energy splittings, AEi w Ai?2, the most 
important (RWA) part of the interaction is 



rr -Aclng) 

-Hint,RWA = -. COS r/i COS ri2PzlPz2 



Ac(ng) . . Ai(ng) 
■ sm rji sm 1J2 



[P+lP-2 + P-lP+2\ 



4 - - 4 

(21) 

This interaction term can be switched off completely only at the symmetry point of at least 
one of the qubits, i.e., when cos 771 = or cos 772 = , by choosing Ug such that the second term 
vanishes. This means, the capacitive coupling constant has to be tuned to a value opposite in 
sign but of the same order as the inductive coupling constant (instead of zero, as one might 
have guessed intuitively). The cancelation is in general not possible for Ej ^ Ec, when Ac (rig) 
is exponentially small. We further note that at the double symmetry point cos 771 — cos 772 = 
the decoupling is stable in linear order with respect to fluctuations of rjk- 

Summary and discussion. - We have analyzed the tunable coupling of qubits mediated 
by a fast coupling element and evaluated the lowest-order nonadiabatic corrections beyond 
the Born-Oppenheimer approximation. When the adiabatic coupling is tuned to cancel an 
additional fixed coupling these nonadiabatic contributions become relevant and can even dom- 
inate. 
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As an example we investigated the capacitive tunable coupler proposed by Averin and 
Bruder, and found that it works only in the regime Ec > Ej. While the total coupling can 
never be switched off completely, this is possible (in RWA) for the most important part of the 
coupling at the degeneracy point of both qubits. In the opposite limit Ej Ec we recovered 
results of an earlier coupling scheme [22] . Here the capacitive coupling (nearly) vanishes, while 
the inductive one dominates. It can not be tuned by the coupler, unless one further modifies 
the design, e.g. in the way as suggested in ref. [7]. 

One strategy to keep the nonadiabatic terms weak is to design a coupler with small adi- 
abaticity parameter |i?qubits|/(-E'i ~ -^'o)) the maximum energy difference between the 
eigenstates of i/qubits should be much smaller than the minimum excitation energy of the cou- 
pler. Another strategy is, to develop designs where [B, iJqubits] = and the Born-Oppenheimer 
approximation becomes exact and, hence, no nonadiabatic corrections appear. Such designs 
have been proposed [8, 13], with tunability achieved via the control of circulating currents in 
SQUID loops. However, the coupling can be switched off only if the system is fabricated with 
identical Josephson junctions in the SQUID loops of the qubits, which is difficult to realize 
experimentally. 
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